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ADDENDUM 

Propagation in a Transport System with 60 Degree Phase Advance. 


The simulations described in this report were performed on an alternate gradient focusing system 
with a low current phase advance of 90 degrees per cell, in this case, the current intensity appears to 
be limited by space charge driven emittance growth. More recent results have shown that a transport 
system with a 60 degree phase advance behaves differently and the current does not appear to be simi¬ 
larly limited. 

By considering the linear stability of a K-V distribution, Laslett and Smith [1] have determined 
that a 60 degree system can stably transport a more intense beam than a 90 degree system. Simulations 
similar to those described here were performed to verify that a K-V system, with £ 60 degree phase 
advance depressed to 24 degrees by space charge, will propagate with no emittance growth or evidence 
of instability. 

More recent simulations performed by Hofmann 12,3] on a 60 degree transport system found that 
a waterbag distribution with a depressed phase advance of 12 degrees showed no substantial rms emit¬ 
tance growth. Motivated by these results, simulations of a K-V distribution with depressed tunes of 12 
degrees and 6 degrees were run. While the K-V distribution does show evidence of the theoretically 
predicted instabilities, the instability saturates at a very low level and rapidly approaches the squared 
distribution described by Hofmann. In contrast to the 90 degree system, very little rapid growth in the 
rms emittance is observed. Some very small slow growth in the rms emittance is observed, but after 
propagation through 100 magnet pairs the rms emittance has grown only about 5%. This absence of 
growth in the rms emittance is similar to what is observed in the simulation of continuously focused 
channels [3]. 

In summary, recent results have indicated that beam propagation in a 60 degree channel appears 
possible at much higher intensities than can be achieved in a 90 degree system. 
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NUMERICAL STUDY OF SPACE CHARGE LIMITED TRANSPORT IN A 
QUADRUPOLE TRANSPORT SYSTEM WITH 90 DEGREE PHASE ADVANCE 

I. INTRODUCTION 

The possibility of using a beam of energetic heavy ions to ignite an inertially confined thermonu¬ 
clear pellet, [1,2] has increased the interest in transporting low emittance high current beams. It has 
been known for some time [3] that such beams are subject to space charge driven instabilities. Ques¬ 
tions concerning the onset of these instabilities and their effect on beam transport have been the sub¬ 
ject of a considerable body [3-9] of theoretical and computational study. 

In this report, several computer simulations are presented which examine the influence of space 
charge driven instabilities on the behavior of a thin lens, alternate gradient, transport system with a 90 
degree phase advance per cell. These simulations follow the evolution of several thousand simulation 
particles in their self consistent electric fields. By following the detailed non-linear evolution of this sys¬ 
tem at various current levels, it is shown that these space charge driven instabilities can cause emittance 
growth. The behavior is found to depend on the current. But the results, which show that increasing 
the current also increases the emittance growth, suggest that this emittance growth, by diluting the 
phase space density of the beam, may place an upper limit on the current intensity which can be tran¬ 
sported. 

In Section II we discuss the simulation results, and in Section III we discuss the implications to 
the transport of intense currents. 

II. SIMULATION RESULTS 

The behavior of quadrupole transport systems in the presence of space charge has been investi¬ 
gated for beams with Kapchinskjj-Vladimirskij (K-V) distributions [11] by integrating the K-V envelope 
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equations 112,13], Stability of K-V systems has in turn been investigated by a perturbation analysis of 
these envelope equation solutions [13-17). Regions of rapid instability, with growth times of just 

several magnet traversals, have been discovered. 

The envelope equation formalism has not, however, yielded insight into the saturation levels or 
the general non-linear behavior of these instabilities. When the amplitude of the unstable modes 
becomes large, the beam can no longer be approximated by a K-V distribution. Though envelope equa¬ 
tions can be used under some circumstances to predict the behavior of non K-V systems [18,19], and 
some recent work has been done on instabilities in non K-V systems [20-22], even the linear behavior 
of these unstable systems is not yet understood. Previous simulations [23] have shown that for 

moderate tune depressions, the system evolves to a steady state after only a factor of two growth in j l 

\ 

emittance. The simulations presented here were performed in order to extend that work to a range of j 

current intensities so as to systematically investigate the non-linear evolution of an unstable K-V sys¬ 
tem. 


The numerical model is a two dimensional particle simulation code of the particle-in-cell type [24] 
optimized for a Texas Instruments ASC computer. The simulation is performed in a plane perpendicu¬ 
lar to the direction of beam propagation, in a reference frame moving with the beam. All variation in 
the beam direction is neglected in this frame. Since the transverse behavior is nonrelativistic, velocities 
and momenta will be referred to interchangeably for notational convenience. 

Numerical tests have been performed [2S] which indicated the existence of a broad numerical 
regime where the simulation results are independent of numerical parameters. These tests showed that 
variations of such parameters ?s time step, grid size, boundary location, and number of particles by fac¬ 
tors of two from the values used, left the results essentially unchanged. The simulations described here 
were generally performed on a 128 x 128 system employing 16384 simulation particles and using six 
time steps to follow the evolution from the center of one magnet to the center of the next. In order to 
eliminate the possiblity that any emittance growth observed is due to numerical elfects, a series of 
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simulations was also performed to test the dependence of emittance growth on the number of particles 
used in the simulation (26] Care was used in all cases to ensure that the number of particles in the 
simulations was sufficient to preclude numerical emittance growth. A thin lens alternate gradient 
transport system with a phase advance of 90 degrees per magnet pair has been assumed for all the 
simulations reported here 

When a system of linear lenses is used to propagate a beam with a K-V distribution, the space 
charge forces are linear and oppose the net focusing forces. All the particles then execute orbits with 
the same phase advance per cell and that phase advance is reduced from the low current value. This 
reduced phase advance is a measure of the amount of current being transported. In fact, for a matched 
beam with a K-V distribution being transported in a thin lens alternate gradient transport system, only 
two parameters are required to completely specify the system. So, for example, specifying the phase 
advance in the absence of space charge and also the reduced or depressed phase advance completely 
specifies the system. Any other way of specifying the total current can also serve instead of the 
depressed tune as a second parameter. 


For a K-V system being transported by finite lenses a measure of the transported current can be 
derived (12] by writing the envelope equations in normalized form. For a thin lens system, however, 
the appropriate [27] measure of beam intensity is: 


Q - 


4 q 1 N r„ f 


(1) 


A /3 y : t 

where q and A are the ion charge state and mass number, N is the number of particles per unit length, 


r p is the classical proton radius, 77 6 is the emittance, and /is the thin-lens focal length. 


Several runs have been performed and these are summarized in Table 1. The cases run are 
identified by the depressed tune, <r. Also tabulated are the values of 0 as defined in Eq. 1. The last 
column in Table I is the number of filaments. A/, formed in the phase space distribution by the 
unstable fields. Note that the first two cases are stable. 
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Figure 1 is a plot of 2048 sample particles (of the 16K in the run) showing four projections of the 
4-dimensional x-Px—y-Py phase space. This plot is of the initial matched K-V distribution of one of 
the runs. A K-V distribution which has the form 



has a uniform projection in each of the four views shown in Figure 1. The graininess of the plots is 
therefore due to the statistical variations in the sample. 

Figure 2 is a similar plot of the phase space of the 45 degree case, but shown after the initial 
equilibrium has passed through 29 magnet pairs. Three well defined filaments can be seen in both the 
x—Px and y—Py projections. Phase space filamentation of this type, which can also be caused by lens 
non-linearities, is known (281 to cause emittance growth. This is also the case here. It is worth noting, 
however, the distinction between non-linearities due to external forces and those seen here. In the 
cases shown here, the non-linearities grow in proportion to their own amplitudes (exponential growth), 
at least in the early stages of the instabilities. Also since the theory of these instabilities, which exam¬ 
ines small perturbations about an equilibrium, does not predict any emittance growth, it is expected that 
any emittance growth will only occur when the unstable perturbation is no longer small. For a particular 
mode therefore, the observation of filamentation and emittance growth is likely to occur near the point 
where the growth of that unstable mode is saturating. 

Figure 3 is a plot of the rms emittances x rm Px rm$ and y rim Py,ms (the cross terms in the normal 
definition [18] are negligibly small) as the the beam evolves from its initial K-V equilibrium. These 
emittances initially remain constant as the unstable modes grow from their small initial values which 
are non-zero due to the statistical imperfections in approximating the initial distribution by a small 
number (in this case 16384) particles. 

Once the initial perturbation has grown to a finite value, the beam phase space density rapidly 
filaments and the emittance grows until a saturated state is reached. The growth then stops. The rapid 
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growth of this mode is in agreement with the linear theory of the stability of K-V distributions in alter¬ 
nate gradient systems. Because only one unstable mode is present for a system with a depressed tune 
of 45 degrees, it has been possible to obtain detailed agreement with the theoretical predictions as to 
growth rate and eigenfunction characteristics for this case |5i. As pointed out previously this agree¬ 
ment occurs only before the major part of the emittance growth is observed 

When the initial current is increased the structure of the most visible unstable mode changes. 
Figure 4 shows the phase space after 20 magnet pairs of an initial K-V distribution with enough current 
to depress the phase advance from 90 degrees to 30 degrees In this case, four filaments are visible. 
Figure 5 shows a continuation of this trend when the tune is depressed further to 15 degrees. Finally, 
Fig. 6 shows the case with tune depressed to 7.5 degrees In this case the filaments are not quite as 
well defined but about 8 to 10 of them are present The more current therefore, the greater the 
number of filaments which are visible. However, the mode which causes the most phase space struc¬ 
ture, and is therefore the most visible, is not necessarily the only one that is unstable. In the higher 
current regime it is in fact likely that several modes are simultaneously unstable. 

Figure 7 is a plot of the v and y rms omittances for the simulation of the case with a 30 degree 

depressed tune. The general behavior, i.e . a period of constancy followed hy a period of rapid growth 

and then a steady state after saturation of this growth is similar to the 45 degree case. The growth rate 

in this case is faster. There is also an amstropy in the emittance in the final state which appears to be 
due to the statistical excitation of the y-growing eigenfunction at a higher amplitude because of the par¬ 
ticular set of random initial conditions on the particles. 

Figure 8 illustrates the occurrence of a transition in behavior when the tune is depressed further. 

The emittance is again a constant initially and then grows rapidly. However, instead of reaching a 

steady state, the emittance continues to grow slowly. Since the distribution function by this lime has 

evolved far from the initial K-V equilibrium, it is difficult to say whether this growth is caused by 

unstable modes which were unstable initially or modes which appeared as a result of the change in the 

distribution function, or in fact, whether the growth is caused by the non-linearities of the distribution 

which resulted after instability saturation However to rule out numerical effects as a cause for this 

5 












I IIABKR 


emittance growlh, several numerical tests were run in which quantities such as the number of particles 
and the spatial resolution of the numerical grid were changed. In none of these cases was the growth 
rate of the emittance affected substantially. 

The simulation run with a depressed tune of 7.5 degrees, as shown in Fig. 9, exhibits an even 
more pronounced persistence to the emittance growth. In order to completely rule out 'he possbility 
that this emittance growth was caused by particle-particle collisions, this case was run using o4o00 parti¬ 
cles. In this case the emittance has grown by about a factor of 1.5 between the end of the rapid growth 
period and the end of the run. This greater emittance growth as the current increases seems to indicate 
the possibility that there is a maximum current which can be stably transported. However, the magni¬ 
tude of this current depends on the details of the scaling of a beam transport system and will be dis¬ 
cussed in the next section. 

The rms emittance is the product of a configuration space and velocity space component, each of 
which can evolve separately. Figure 10 shows, for the case with 45 degree depressed tune, a plot of the 
x-extent of the beam at the center of an x focusing lens as the system evolves. Three curves labelled r, 
a. and m at the right margin correspond to the rms, absolute average and maximum (the position of the 
maximally displaced particle) measures of the size. 

Figures 11-13 show the same graphs for the runs with depressed tunes of 30, 15. and 7.5 deg r ees 
respectively Both the rms and absolute average beam sizes change little. Therefore little of the growth 
in the rms emittance shown in Figures 3, 7, 8 and 9 occurs due to growth in the configuration space 
component. In all cases the maximum particle, or edge of the distribution function, has moved to 
about twice its initial value. The beam size has not increased much for most of the particles, but some 
particles have been scattered out to twice the initial edge of the beam. 

Figures 14-17 are the corresponding graphs of the evolution of the beam extent in P.v space, also 
taken at the center of an x-focusing lens. These curves also illustrate that most of the growth, espe¬ 
cially at the highest currents, is due to growth in the velocity component of the emittance. The edge of 
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the beam has moved way out in velocity space. Since the same spreading has not occurred in 
configuration space, this is an indication of the highly non-linear space charge forces near the edge of 
the beam, with the non-linearity becoming more pronounced at the higher currents. 

These curves all refer to the behavior at the center of the x-focusing lenses. This was done 
because this is where the beam is at its maximum in the x direction, and it is therefore appropriate to 
look at the beam growth at this point in order to determine the size channel which would be necessary 
to propagate the beam. 

The self-consistent equilibria corresponding to a K-V distribution are not difficult to visualize. All 
projections of the beam are uniform ellipses. The non K-V steady states which the system has evolved 
to in these simulations are somewhat more complex. Figure 18 shows four views of the x-Px—y-Py 
phase space from the simulation with an initial 30 degree depressed tune after the beam has propagated 
through 100 magnet pairs. Figures 19 and 20 are contour plots of the x-Px and y-Py phase spaces of 
Figure 18. Adjacent contours are separated by a factor of four in density. This run was chosen because 
it was the highest current at which a steady state was reached. The shapes of the x-Px and y-Py con¬ 
tours differ because the plot was done in the center of an x-focusing lens. A similar contour plot done 
at the center of the y-focusing lens shows the expected reversal of the two plots. 

Another set of plots can be generated of the configuration space and velocity space contours. 
These are most conveniently taken half way between the x-focusing and t-focusing lenses. At this 
plane a matched equilibrium is symmetric in the x and y directions. Figures 21 and 22 show these plots 
for configuration and velocity space respectively for the same case run with a 30 degree depressed tune 
taken at the end of the simulation. 

The squaring of these contour lines is just one indication as to the difficulty of analytic characteri¬ 
zation of the high space charge equilibria. This squaring increases as the current is increased. Since the 
contours in phase space are not ellipses, it is also an indication that any definition of emittance in terms 
of the percentage of particles within an elliptical phase space contour may not be appropriate. 
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For completeness, it is also possible to plot various cross sections of these curves. Figure 23 is 
such a cross section of Figure 21. Figure 23 is a plot of fix), the number density as a function of x, at 
the y = 0 plane integrated over all velocities. Because this curve is generated at the symmetry plane, it 

is the same form as a curve of f(y) at x - 0. Figure 24 is a similar plot of f(Py) at Px = 0 integrated 
over all positions. These curves are typical. A reasonable approximation to them can be obtained as 

the sum of a distribution with a parabolic profile and a smaller one with a linearly decreasing tail (21] 
The number of such curves which can be generated for all the cases run is large. Most of the informa¬ 
tion contained in these curves can be summarized by noting that as the current increases, the contours 
get squared as do the configuration space cross sections. The velocity space cross sections develop 
increasingly long tails. The information of the length of these tails however, can be pretty satisfactorily 
estimated from the curves of the evolution of velocity vs. time (Figs. 14-17), by comparing the max¬ 
imum velocity particle to the rms velocities. However at the mid plane the velocity and configuration 
behavior look generally similar to Figs. 23 and 24. 

III. CONCLUSIONS 


Current intensity depends on the ratio of current to emittance. Therefore as the current 
increases, if the emittance increases correspondingly, then the depressed tune is unchanged. In fact, if 
the distribution function remains unchanged, the beam is then represented by the same solution to 
Maxwell’s equations in the transverse direction. It is possible from purely dimensional arguments to 
show that any solution to Maxwell’s equations in the transverse direction corresponds to a familv of 
solutions at different currents. These currents can be related to the power propagating in a transport 
system by the relation 


- C — 
z 


4/3 

(B q t) i,} (fiy) Vi (y-\) (3) 

where P is the power in watts, A and z are the mass number and charge stale. B q is the quadrupole pole 
tip field in Tessla, and rr « is the emittance in meter-radians. The coefficient C contains the informa¬ 
tion on the current intensity and the beam distribution function. 
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This formula was first suggested by Maschke and was derived by Courant 129] for a K-V distribu¬ 
tion. A different treatment was derived by Reiser [30]. The derivation in terms of purely dimensional 
considerations is presented in the appendix. 

It is possible to determine a value for the coefficient C in Eq. (3) corresponding to any solution to 
Maxwell’s equations in general, and therefore specifically for the distributions coresponding to each of 
the simulations performed. However, the power transported is expressed in Eq. (3) in terms of the 
pole tip field. Therefore because the simulations were performed using thin lenses, and these are 
unphysical as they require singular pole tip strengths, some assumption must be made about the magnet 
fill in the transport system. From the extensive linear stability calculations [13-1S], it is known that, at 
least in the linear stages of the instability, the behavior of the system is largely insensitive to magnet fill 
as long as the phase advance per cell is kept unchanged. This insensitivity to lens fill has also been 
seen in simulations [7], 

Laslett has suggested [31] the definition of a beam figure of merit, / = (?/<w max ) 2/ \ based on the 
behavior of the scaled envelope equation (121. The coefficient C in Eq. (3) is related to F by, C = 3.43 
x 10 15 fin MKSA units. For a K-V system, the values of F have been calculated for a range of param¬ 
eters [31-33]. These results are plotted in Fig. 25, as a function of depressed tune, for a system with 
90 degree phase advance and with 50 percent lens fill. The values of / can be calculated beyond the 
range shown in Fig. 25 from the expression [31] F = 9.5 o- _2/3 , which is accurate to within one percent 
below 10 degree depressed phase advance. From the K-V values, it is possible to calculate the value of 
F, and therefore the coefficient C, appropriate to the final state in each of the simulation runs. 

Since the current does not change during the simulation, the coefficient in Eq. (4) changes due to 
the growth in emittance as the beam evolves. Because the beam size also increases an additional factor 
must be included to account for the fact that the magnet pole lip must be moved out. These values are 
listed in Table II. The value at 60 degrees, which is not unstable, is the same as Laslett's value. The 
other values are calculated from Laslett's values by dividing by the increase in rms emittance to the 2/3 
power and then also dividing by the increase in size to the beam edge (as measured from Fig. 13-16) 
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also raised to the 2/3 power. This second factor is necessary if the pole tip is moved out far enough to 
include the entire beam. Since the emittances used are rms emittance these particular coefficients are 
valid only when using this rms measure. 

The last two entries in the table correspond to distributions which are not in steady state, but are 
still evolving. Notice for example that the initial figure of merit increases by a factor of 1.6 between 15 
and 7.5 degrees, by the end of the run this improvement has been reduced to 1,2. If the run were con¬ 
tinued further, this reduction would have continued still further resulting in the possibility of no gain in 
transported current intensity. 

For the simulations performed, that is a constant current propagating down a uniform transport 
line, a K-V system can become rapidly unstable. This rapid instability appears to saturate after about a 
factor of two growth in emittance. At the highest current intensities studied, this period of rapid 
growth is followed by a period of slower growth which tends to limit the current intensity than can be 
propagated over long distances. 

Because of the inherently non-linear nature of the phenomena involved it is difficult to say any¬ 
thing about beam distributions more general than K-V If it as assumed that the behavior of non K-V 
distributions is adequately represented by the behavior of the simulations alter initial saturation, then it 
would be possible to propagate very intense currents for modest distances without too much emittance 
growth However the existence of modes which can grow rapidly and which could be excited if the 
beam distribution would evolve to a state closer to a K-V distribution, suggests that caution be exer¬ 
cised. Conservative design still seems to preclude operation in this space charge limited regime without 
a much fuller understanding of the effect of beam distribution function on stability. 
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Table I. Summary of Cases Run 


cr 

0 

M 

90 

0.0 

— 

60 

0.93 

— 

45 

1.6 

3 

30 

3.1 

4 

15 

6.2 

6 

7.5 

12.4 

8-10 


Table II. Transport Figures of Merii Modified 
by Space Charge Driven Instabilities 


<T 

F, 

Ft 

60° 

0.34 

0.34 

45° 

0 56 

0 22 

30° 

0.87 

0.31 

15° 

1 52 

0 54 

7.5° 

2 48 

0.67 


IS 




Appendix 

CURRENT SCALING IN A STRONG FOCUSING CHANNEL 


Consider a beam of charged particles, in the paraxial approximation, in a reference frame moving 
with the beam. The transverse dynamics, in this frame, is described by the Vlasov-Poisson system of 
equations. 

+ v-V/ + (E et + E) • V v / = 0 (Al) 

where the distribution function is normalized so that 

J* / dv x dv v dx dy = 1. 

and E,. v represents the focusing magnet force in the beam reference frame. The space charge fields are 
found from 


V • E = 4fl-p. (A2) 

In order to investigate the scaling of this system we write these equations in terms of normalized 
variables. Setting r = r 0 r, v = v„ v and t = t„i with r„ = v„/„, and also E,,, = E 0 E ex , 

H + v • V/+ ~~ (£„ E (><( + E) • V, / = 0. ( A3) 

dr v,fm 

Since the divergence operator is just a linear operator for the purpose of this scaling analyses, it 
can be formally inverted and written as: 


E = 4 rrLp (A4) 

where L is just the appropriate linear operator inverse to the divergence. The self-consistent electric 

field can then be written in terms of the total current / as 


E = 


47T / r„ L 
r oPyc 


f f d\\ d\\ . 


(A5) 
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Equation (A3) may then be rewritten as: 


H + ‘"V/ + (C, L, + C 2 L / / </v v </v,) ■ V, / = 0 


where 


- -T 2 - £„ 


P ycv„m 


Any distribution function /satisfy (A6) with a particular form of £„, will therefore be a solution 
for any combination of parameters which leave C\ and C 2 unchanged. We can use this to derive a scal¬ 
ing relation for /. 


Introducing an emittance 


e = C, r 0 vjfiyc. 


mC i v,,m 


Solving for the current 


— e 2/ 3 E** y V 3 


Cf’ 4v f 


If the external focusing force is written in terms of transverse magnetic fields, E„ = PyB 0 , this 


becomes 


/ = C\,j^ ' f 2/} B* 3 (py) Vi . 


The power transported by this current, is 


P = (y-1) — / 
<1 


e 2/ ' fl/ 3 (0y) 7/3 (y-l) 
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0 x -•< py t. 

Fig ' — Four projections of the four-dimensional K-V distribution showing 2048 sample particles. The configuration 
space units are numerical grid cell and the velocity (momentum) units are in grid cells traversed in one magnet period. 


18 







p 

X 




l*«« 


ttl -|4 


P Y 


14 


Four phase space projections showing filament formation in phase space after traversing 29 magnet pairs The 
initial dish ibution is a K-V system with a depressed phase advance of 45 degrees. 
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I ig ' - Evolution of the x and v rms emlttances in a system with an initial matched K-V distribution with enough 
current to depress the tune from 90 degrees to 45 degrees. The vertical scale is normalized to the initial emittance. 
The horizontal scale is measured in alternate gradient magnet pairs. 








f ig -I — Four pha.se spate projections showing the evolution, after 20 magnet pairs, of a K-V system 
with a depressed phase advance of 20 degrees 
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Evolution of the x and y rms emiltances in a system with an initial K-V distribution with a depressed phase 

advance of 30 degrees 




Evolution of the x and y rms emittances in a system with an initial depressed phase advance of 15 degrees. 
Some slow long time scale growth in the emittance is apparent. 









Evolution of the x and > rms emittances in a system with an initial depressed 
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Evolution of the beam size in the x direction for an initial 15 degree phase advance 










simulation vuth a 4^ degree initial phase advance 
















( volution of the v-momenta for the vase with ^ 5 degree phase advance 









Ill * *< 


I tg in — Four views of the x-Px-y—Py phase space of a simulation with an initial depressed phase advance of 30 de¬ 
grees which has been propagated through 100 magnet pairs. The system has gone unstable and by this point has 
reached a steady state 
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lit i 1 - A contour plot of the v- l\ phase space projection shown in fig IS Adjacent countours are separated bv a 

factor of four in density 
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A contour plot of the density in real space of the simulation with 30 degree depressed phase advance after 
100 magnet pairs The plot is taken midway between an v foeuslng and r-focusing lens 
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